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ABSTRACT BASIC SPHERICAL CLUSTERS 

Packed spherical particle beds have wide applications throughout Hexagonal or cubic densely packed beds have 12 neighbors and a 

the process industry and are usually analyzed using an appropriate void fraction of about 0.26, Fig. L Two basic spaces can be identified: 
combination of laminar and turbulent flows with empirically derived 

coefficients of which the Ergun (1952) relation is probably the best Tetra-space. — The unit is formed of four (4) spheres with centers 

known. The 3-D complexity of the void distribution within the beds at the apex of a tetrahedron. Three spheres lie in a plane touching one 

has precluded detailed studies of sphere clusters. Numerical modeling another with the centers forming an equilateral triangle. The fourth 

and flow vector visualization for the basic tetra- and hexa-sphere sphere is placed in the hollow formed by the three spheres forming a 

clusters and two hyper-sphere clusters are presented at two Reynolds regular tetrahedron. Thus there are four (4) open curvilinear triangular 

numbers, 400 and 1 200. Cutting planes are used to enable visualization faces for flow and four (4) spherical triangular surfaces or boundaries, 

of the complex flows generated within the sphere clusters and are Fig. 2(a). Alternately, a half tetra-space unit is the three spheres in a 

discussed herein. The boundary conditions and flow fields for the plane with a sphere placed in the hollow on each side and the plane of 

simple clusters are also compared to the hyper-clusters with larger symmetry passing through the centers forming an equilateral triangle, 

variations found for hexa- clusters. Hexa-space. — The unit is formed from six (6) spheres. Three 

spheres lie in one plane touching one another with centers forming a 
equilateral triangle and three others lie in a second plane stacked upon 
INTRODUCTION the first plane and rotated 60° to fit the hollow of those in the first plane 

emulating a six pointed star. An alternate method of constructing this 
Although flow in packed beds of spheres may be correlated in the space is to arrange four spheres in a plane at corners of a square, and 

same manner as flows through brush seals and arrays of pins (cylin- place two spheres in the hollows of the layer one on each side. The 

ders) via the simple combination of laminar and turbulent flows with centers of the spheres form a double pyramid. To visualize this, rotate 

appropriate coefficients (e.g., Bird et al. (1960), Hendricks et al. the hexa-space, figure 1, about the solid horizontal line (upper layer-2), 

(1996)) yet flow details have not been delineated. The complexity of away from the viewer, until the solid and dashed horizontal line (lower 

any useful bed of spheres includes occlusions, voids, irregular and layer- 1) coincide. In this view the double pyramidal structure becomes 

different diameter spheres, compaction strains, crushed and cracked apparant. This space has eight (8) open faces for flow and six (6) 

spheres, debris holdup and variations in thermophysical properties. At spherical rectangular surfaces or boundaries, Fig. 2(b). 

this point, these factors will be considered secondary in an attempt to 

calculate and visualize the flow in a regular bed. More details on flows Subdomains of the tetra and hexa Interstitial Spaces 
in regular beds of spheres can be found in Athavale et al. (1994a). 

Visualization of flows in packed beds of spheres as described herein To generate subdomains that can be repeated, a common shared 

is comprised of(l) mathematical modeling and numerical representa- point(centroid) edges and surfaces are required. The centroid of 

tion of the flow fields and (2) physical modeling and laser sheet interstitial space was taken as this point, with straight lines between the 

illumination of flow field particles. Unlike visualization of flows in centroid and the corners of the open faces defining the shared edges 

pinned arrays and brush seals which are weakly 3-D dependent, flow, and shapes of the subdomains. These subdomains are either tetrahe- 

heat and mass transfer in spherical beds are strongly 3-D dependent dral or pyramidal in shape, and are defined either with the open area 

and require the combined efforts of numerical and experimental or the curved surfaces of the sphere as the basis that is joined to the apex 

techniques to embody and characterize models and results. located at the central point. 
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Tetrahedral- 1. — This subdomain has the triangular curved sphere 
surface as the base Fig. 3(a). The base shape is generated in tetra-space. 
The corners of the triangle are joined to the apex to form the subdomain 
and is needed only for tetra-spaces. 

TetrahedraI-2. — This subdomain has the flat, open area shown in 
Fig. 3(b) as the base, and the three corners are joined to the central apex 
to form the subdomain (Fig. 3(c)). Since the subdomain has the open 
cutting plane as the base it is required in both tetra- and hexa-spaces 
when forming hyperclusters. 

Pyramidal. — This subdomain has a quadrilateral portion of the sphere 
surface as the base, and the four comers are joined to the apex to form 
the subdomain, Fig. 3(d) and required only in the hexa-space. 

A compilation of the subdomains and the corresponding surface grids 
are shown in Figs. 3(e) to 3(g). The axis convention for grid generation 
are illustrated in Fig. 4(a) for tetrahedral- 1, Fig. 4(b) for tetrahedral-2 and 
Fig. 4(c) for pyramidal. 

The geometries of the spaces are complex, and a global axis conven- 
tion as well as numbering sphere locations was found to be necessary and 
convenient in describing the geometry as well as analyzing the flow- 
field. This convention is shown in Figs. 5(a) and (b). 

The interstitial spaces are created using a group of the subdomain with 
the apexes coincident at the space centroid. Tetrahedral domains are used 
for all open faces and walls in the tetra-space with pyramidal domains 
used for walls in the hexa-space. All apexes coincide at the center of the 
interstitial space which implies a polar point. Tetrahedral domains also 
have a polar line or centerline. Each tetra-space has 8 domains and each 
hexa-space has 14 domains. 

The tetra-space has tetrahedral- 1 subdomains on the sphere surfaces 
and four tetrahedral-2 subdomains on the four open faces, Fig. 6(a). 
The hexa-space had six surrounding spheres, with bases on the sphere 
surfaces. Eight cutting planes with eight open area used as bases for the 
tetrahedral-2 subdomains, Fig. 6(b). 

The packed bed is generated by populating a region with the tetra- and 
hexa-spaces joined at the open faces in each. Care is needed to ensure 
that the resulting clusters generate the original sphere surfaces correctly. 
Two types of first-level hyperclusters are possible: 

Hypercluster- 1. — Tetra-space connected to 3 hexa-spaces and one 
plane of symmetry with 54 domains, Fig. 7(a). 

Hypercluster-2. — Flexa-space connected to 8 tetra-spaces with 78 
domains Fig. 7(b). 

An algebraic grid generator was used to generate clustered grids in 
each of the subdomains which were then assembled to form multiblock 
grids for the interstitial spaces and hyperclusters. The tetra-space had 
30(i)xl0(j)xl0(k) for each domain (24000 cells), while the hexa-space 
had 8x8x8 pyramid domain, 24(i)x8(j)x8(k) tetrahedral domain 
(15 360 cells). The clusters had 24x8x8 for the tetrahedral domain and 
8x8x8 for the pyramidal domain. 

Flow Solver 

The numerical simulations were performed using SCISEAL, a 
pressure based, finite volume 3-D CFD code with implicit multiblock 
capability, Athavale et al. (1994b). The working fluid was nitrogen gas 
at pressures of 1.5 MPa. The flow Reynolds numbers were 400 and 
1200 and the flow was assumed to be laminar with constant gas 
properties. The boundary conditions for the single tetra- and hexa- 


spaces are shown in Figs. 6(a) and (b). Similar conditions were also 
used for the two hyperclusters. 

RESULTS 

Tetra-space 

Included here are several representative velocity vector plots 
for two Reynolds numbers 400 and 1200, taken along selected 
cutting planes. Figure 8(a) shows the velocity vectors at x - 1.1 cut- 
ting plane, which is a little to the right of the symmetry plane at x = 1 .0, 
see Fig. 5(a). The uniform, parallel inlet velocity is seen on the left 
edge of both plots in Fig. 8(a). The frontal stagnation zone created by 
the incoming flow on sphere 3 is seen on the right side, and the large 
recirculation region at the bottom of the figure is created by a 
combination of the stagnation on sphere 3 and the symmetry boundary 
on the bottom plane of the tetra-space. As expected, this recirculation 
zone is somewhat stronger in the high Reynolds number flow. In the 
upper portion, the flow splits and moves upwards (in +z direction) and 
sideways (± x-directions) as it exits through the open areas on either 
side of sphere 3. There is evidence of separation of the flow on the 
upper wall, but the strength of this recirculation zone is nearly the same 
for both Reynolds numbers. 

The incoming flow generates separation bubbles near the surfaces 
of spheres 1 and 2 (due to the surface curvature); to see these zones, 
vector plots along the z = 0.18 cutting plane are shown in Fig. 8(b). 
The largest contrast between the two Reynolds number flows is seen 
in this figure. A double vortex is seen in the plot for the higher 
Reynolds number flow, which is nearly absent in the low Reynolds 
number case. The effect of this bubble would have a direct impact on 
surface temperature distributions for example. The frontal stagnation 
zone on the surface of sphere 3 is seen in both the plots. 

Vector plots along the y = 0.63 cutting plane are shown in 
Fig. 8(c) and part of the separation zones on sphere 1 and 2 are seen. 
The incoming flow is out of the plane of the paper. The flow near the 
sphere surfaces in both cases turns upwards as it moves toward the exit 
boundaries on either side of the plots. The recirculation bubbles for the 
higher Reynolds number flows are again seen to be a little stronger. 

Hexa-space 

The hexa-space flow is more complex due to multiple inlet 
boundaries, Fig. 6(b). The velocites at the x = 1.0 cutting plane (a 
symmetry plane) are shown in Fig. 9(a). The front stagnation region 
on sphere 3 is seen in both plots and a very weak recirculation in the 
central portion. A recirculation zone near the open bottom surface is 
seen, and is nearly the same for both Reynolds numbers. This 
recirculation is also seen to draw fluid in through the bottom face. 
Flow near sphere 3 turns upwards as it exits though the open face on 
the right. A central vortex is seen on the symmetry boundary in both 
cases, and is a result of the upward flow motion on the right side, 
which forces the flow to go in the negative z direction under the 
surface of sphere 4. The incoming flow on either side of sphere 4 
produces a double vortex under this sphere, which is seen in the 
z = 1.04 section, Fig. 9(b). The double vortex is somewhat stronger 
and smaller in the higher Reynolds number case. This double vortex 
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forces the flow in the negative y and z directions under sphere 4, and 
this was also seen in the previous figure. At a lower section z = 0.36, 
Fig. 9(c), the incoming flow generates symmetric vortices near the 
surfaces of spheres 1 and 2 as in the case of the tetra-space, and the 
vortices are stronger for the higher Reynolds number. Vector plots at 
y = 0.58 cutting plane are shown in Fig. 9(d). The double vortex is 
seen under sphere 4 in both plots. The flow is out of the plane of the 
paper. A small recirculation is seen near the top open face which 
indicates some inflow as well as some outflow in this face. The 
overall structure of the flow, except the central double vortex is quite 
similar for the two Reynolds numbers. 

Hypercluster -1 

The central cutting plane at x = 1.01 Fig. 10(a) shows the incoming 
flow at the left, which now show some variation in the velocity 
magnitudes (Reynolds number of 1200), including some evidence of 
small recirculations near the x-y plane, where a sphere surface from the 
adjoining hexa-space is present. The recirculation bubble in the central 
tetra-space near the bottom edge is weaker and smaller in size as a 
result of the inclination of the incoming flow through the plane on the 
left. There is also some evidence of recirculation near the top edge, as 
the flow exits through the right-side boundary. 

The vector plot along, z = 0.18 shown in Fig. 10(b) has some 
differences from the single molecule plot Fig. 8(b). Although the 
velocity pattern near the spheres 1 and 2 is similar, the hyper-cluster 
solution shows an absence of the double vortical structure seen in 
Fig. 8(b.2). The y = 0.63 plane, Fig. 10(c) shows a flow pattern similar 
to that for the single tetra-space Fig. 8(c.2). Several recirculation zones 
exist in the region attached to the central interstitial space shown in 
these figures. The solutions on the single spaces were obtained using 
a specified constant pressure on the exit boundaries. The pressure 
values on the corresponding interfaces in the central tetra-space were 
checked to assess the validity of this assumption. The overall pressure 
variation across these faces was found to be sufficiently small to justify 
the constant value assumption used in the single space calculations. 
Pressure plots were developed but not reported herein. 

Hypercluster-2 

Six pyramidal subdomains in the central hexa-space and 72 tetrahe- 
dral subdomains for a total of 113664 cells. The velocity vectors 
(Reynolds number of 1200) along the plane at x = 1.0 are shown in 
Fig. 11(a). As in the previous case some differences from its single- 
space counterpart Fig. 9(a.2), are noted. The incoming flow from the 
left edge has a faster portion near the top edge which is a result of the 
tetra-space present on the left of this interface. This jet tends to destroy 
the weak central recirculation zone seen earlier. The recirculation zone 
near the bottom interface is still present, and the double vortical 
structure under sphere 4 is seen. The plane along z = 1.04 also shows 
this double vortex clearly in Fig. ll(b.l). However, the vortex is 
somewhat weakened due to the change in the inlet velocity profiles at 


the interface boundaries on either side of sphere 4. The section along 
z = 0.36 in Fig. 1 1 (b.2), however shows a nearly smooth flow along the 
spheres 1 and 2 while the single space counterpart Figs. 9(b.2) and (c.2) 
show a vortex near the sphere surfaces. As in the case of hypercluster- 
1 , this is due to the inclination of the incoming flow. The double vortex 
under sphere 4 is also seen in the y = 0.58 plane shown in Fig. 1 1(c). 
This figure also shows the changes that are present in the inlet flows 
on either side of sphere 4, where the flow exits tetra-spaces attached 
to the central hexa-space, compare with Fig. 9(d.2). 

SUMMARY 

Two basic interstitial flow geometries, tetra- and hexa-spaces, have 
been identified as essential elements for generating hexagonal or cubic 
densely packed ideal beds. The interstitial spaces were divided into 
suitable subdomains for structured grid generation and subsequent 
solutions. The flow solver SCISEAL, a pressure based multidomain 
Navier-Stokes solver was used to obtain representative flow solutions 
on the tetra- and hexa-spaces as well as first-level hyperclusters 
generated using the basic spaces. 

The flow is fully 3-D and complex with multiple vortical structures 
that generate local stagnation, separations, and most importantly flow 
threads or fast flowing regions in the direction of the bulk flow field. 

Preliminary solutions on first level hyperclusters indicate that the 
results differ from the basic structures; still the major vortical struc- 
tures and flow threads remained. 

The vortical structures and local stagnation zones are inferred to 
have a significant influence on transport phenomena. 

The results, although not compared herein, are in good qualitative 
agreement with experimental observation. 

REFERENCES 

Athavale, M.M., Przekwas, A.J., Braun, M.J. (1994a) Flow and Heat 
Transfer Analysis for Packed Beds on Heated Spheres. Part 1: 
Numerical Computations. Final Report Contract NAS3-25967. 
NASA Lewis Research Center, Cleveland Ohio. 

Athavale, M.M., Przekwas, A.J., Hendricks, R.C., and Liang, A., 
(1994b) A SCISEAL: A Three-Dimensional CFD Code for Accu- 
rate Analyses of Fluid Flow and Forces in Seals, 1994 Conference 
on Advanced Earth- to- Orb it Propulsion Technology, NASA 
Marshall Space Flight Center, May 1994. 

Bird, R.B., Stewart, W.E., and Lightfoot, E.N.(1960), Transport 
Phenomena, John Wiley Press, New York, p. 200. 

Ergun, S. (1952) Chem Engr. Prog. 43(92). 

Hendricks, R.C., Kudriavtsev, V.V., Braun, M.J., Athavale, M.M. Flows in 
Pinned Arrays Simulating Brush Seals - to be presented at the hit Cong. 
On Fluid Dynamics, Dec. 29-31, 1996 Cairo Egypt 
Kerrebrock, J.L. and Kalamas, J. (1993) Flow Instability in Particle- 
Bed Nuclear Reactor, AIAA-9 3-1758. 


3 
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Figure 1 .—Stacking of sphere layers to generate densely-packed arrangements. 
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Figure 2.— (a) Solids model for the tetra-space. (b) Solids model for hexa-space. 
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Figure 3.— Outlines of the three basic subdomains used in the tetra- and hexa-spaces. (a) Tetrahedral-1. (b) The 
cutting plane used to define the outer boundaries of the two clusters. View along the normal to the cutting plane, 
(c) Tetrahedral-2. (d) Pyramidal, (e) Grid for tetrahedral-1 subdomain with sphere wall base in tetra-space. (f) Grid 
for tetrahedral-2 subdomains used in the tetra-and hexa-spaces. (g) Grid for pyramidal subdomain used in 
hexa-space. 
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Figure 4(a)— Surface grids for tetrahedral-1 subdomain with axis convention. (1) Grid on the sphere surface. 
(2) Grid on one of the triangular sides. 
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Figure 4(b) — Surface grids for tetrahedral-2 subdomain with axis convention. (1) Grid on the open base area. 
(2) Grid on one of the triangular sides. 
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Figure 4(c)— Surface grids for pyramidal subdomain with axis convention. (1) Grid on the sphere surface. 
(2) Grid on one of the triangular sides. 
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Figure 5.— (a) Axis convention and sphere arrangement for the tetra-space. (b) Axis 
convention and sphere arrangement for the hexa-space. 
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Nomenclature 
[1,4,6] =plane cutting 

sphere 1,4, and 6 
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Figure 6(a). — Solids model, flow geometry and boundary conditions for the single tetra-space. 
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Nomenclature 
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to x-y plane preSSUreS 



K [1,3,6]/ \- 

- [3,5,6] behind 


m/1 

sphere #5 


■Lr 

z\ 


■Ik 

\ 

/Sphere #5 


Sphere #4 
surface grid 
shown for 
clarity -|[ 




Inlet 

flow 



Specified 

pressures 

[2,3,5] 

behind 

spheres 

3 and 2 


, / l2i 

V [1,2 ,4U 2 ' 4 


Uniform inlet 
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Sphere #1 

Figure 6(b). — Solids Model, flow geometry and boundary conditions for the single hexa-space. 
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Figure 7(a).— Perspective view of hypercluster- 1. Tetra-space in center, three hexa-spaces surrounding it. 
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space #1 


Figure 7(b). — Per-perspective view of hypercluster-2. One hexa-space in center, eight tetra-spaces 
surrounding it. 
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Figure 8(a). — Velocity vectors in the single tetra-space along the x = 1 .1 cutting 
plane. (1) Reynolds number = 400. (2) Reynolds number = 1200. 
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